function hp(x)
y1=0.0223*x^3-1.4408*x^2+61.988*x+9.5765
y2=0.0044*x^3-0.2786*x^2+27.0449*x+115.9317
y3=0.0011*x^3-0.0579*x^2+10.5959*x+39.6062